R version 3.1.2 (2014-10-31) -- "Pumpkin Helmet"
Copyright (C) 2014 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin10.8.0 (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> print(formulae)
[[1]]
TP ~ W * L + Error(years)
<environment: 0x1028e0c28>

[[2]]
Chl ~ W * L + Error(years)
<environment: 0x1028d26e8>

[[3]]
Therm ~ W * L + Error(years)
<environment: 0x1028ca2d8>

[[4]]
Epi ~ W * L + Error(years)
<environment: 0x1028bd8e8>

[[5]]
Light ~ W * L + Error(years)
<environment: 0x1028b41c8>

[[6]]
Ice ~ W * L + Error(years)
<environment: 0x102899780>
> print(mlist)
$TP

Call:
aov(formula = formulae[[i]], data = d2015)

Grand Mean: 441.8475

Stratum 1: years

Terms:
                Residuals
Sum of Squares   46650.69
Deg. of Freedom        15

Residual standard error: 55.76778

Stratum 2: Within

Terms:
                       W        L      W:L Residuals
Sum of Squares    872.66  1361.14    21.81  48637.92
Deg. of Freedom        1        1        1        45

Residual standard error: 32.87617
Estimated effects may be unbalanced

$Chl

Call:
aov(formula = formulae[[i]], data = d2015)

Grand Mean: 10.42045

Stratum 1: years

Terms:
                Residuals
Sum of Squares   89.85201
Deg. of Freedom        15

Residual standard error: 2.447475

Stratum 2: Within

Terms:
                        W         L       W:L Residuals
Sum of Squares   0.022775  8.793033  0.064183 21.093732
Deg. of Freedom         1         1         1        45

Residual standard error: 0.6846529
Estimated effects may be unbalanced

$Therm

Call:
aov(formula = formulae[[i]], data = d2015)

Grand Mean: 8.724916

Stratum 1: years

Terms:
                Residuals
Sum of Squares   1.706237
Deg. of Freedom        15

Residual standard error: 0.3372672

Stratum 2: Within

Terms:
                        W         L       W:L Residuals
Sum of Squares  0.0131248 0.0439710 0.0934546 2.1784973
Deg. of Freedom         1         1         1        45

Residual standard error: 0.2200251
Estimated effects may be unbalanced

$Epi

Call:
aov(formula = formulae[[i]], data = d2015)

Grand Mean: 8.636609

Stratum 1: years

Terms:
                Residuals
Sum of Squares   2.874749
Deg. of Freedom        15

Residual standard error: 0.4377784

Stratum 2: Within

Terms:
                       W        L      W:L Residuals
Sum of Squares  0.230754 0.355258 0.028031  3.253536
Deg. of Freedom        1        1        1        45

Residual standard error: 0.2688881
Estimated effects may be unbalanced

$Light

Call:
aov(formula = formulae[[i]], data = d2015)

Grand Mean: 1.217145

Stratum 1: years

Terms:
                 Residuals
Sum of Squares  0.02130495
Deg. of Freedom         15

Residual standard error: 0.03768727

Stratum 2: Within

Terms:
                          W           L         W:L   Residuals
Sum of Squares  0.000122095 0.002320080 0.000036659 0.007045864
Deg. of Freedom           1           1           1          45

Residual standard error: 0.01251298
Estimated effects may be unbalanced

$Ice

Call:
aov(formula = formulae[[i]], data = d2015)

Grand Mean: 0.02466637

Stratum 1: years

Terms:
                  Residuals
Sum of Squares  0.001886981
Deg. of Freedom          15

Residual standard error: 0.011216

Stratum 2: Within

Terms:
                           W            L          W:L    Residuals
Sum of Squares  0.0010158816 0.0000040508 0.0000106957 0.0014317304
Deg. of Freedom            1            1            1           45

Residual standard error: 0.005640588
Estimated effects may be unbalanced

> print(lapply(mlist, summary))
$TP

Error: years
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals 15  46651    3110               

Error: Within
          Df Sum Sq Mean Sq F value Pr(>F)
W          1    873   872.7   0.807  0.374
L          1   1361  1361.1   1.259  0.268
W:L        1     22    21.8   0.020  0.888
Residuals 45  48638  1080.8               

$Chl

Error: years
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals 15  89.85    5.99               

Error: Within
          Df Sum Sq Mean Sq F value  Pr(>F)    
W          1  0.023   0.023   0.049   0.827    
L          1  8.793   8.793  18.758 8.2e-05 ***
W:L        1  0.064   0.064   0.137   0.713    
Residuals 45 21.094   0.469                    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

$Therm

Error: years
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals 15  1.706  0.1138               

Error: Within
          Df Sum Sq Mean Sq F value Pr(>F)
W          1 0.0131 0.01312   0.271  0.605
L          1 0.0440 0.04397   0.908  0.346
W:L        1 0.0935 0.09345   1.930  0.172
Residuals 45 2.1785 0.04841               

$Epi

Error: years
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals 15  2.875  0.1916               

Error: Within
          Df Sum Sq Mean Sq F value Pr(>F)  
W          1  0.231  0.2308   3.192 0.0808 .
L          1  0.355  0.3553   4.914 0.0317 *
W:L        1  0.028  0.0280   0.388 0.5367  
Residuals 45  3.254  0.0723                 
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

$Light

Error: years
          Df  Sum Sq Mean Sq F value Pr(>F)
Residuals 15 0.02131 0.00142               

Error: Within
          Df   Sum Sq   Mean Sq F value   Pr(>F)    
W          1 0.000122 0.0001221   0.780 0.381899    
L          1 0.002320 0.0023201  14.818 0.000371 ***
W:L        1 0.000037 0.0000367   0.234 0.630822    
Residuals 45 0.007046 0.0001566                     
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

$Ice

Error: years
          Df   Sum Sq   Mean Sq F value Pr(>F)
Residuals 15 0.001887 0.0001258               

Error: Within
          Df    Sum Sq   Mean Sq F value   Pr(>F)    
W          1 0.0010159 0.0010159  31.930 1.03e-06 ***
L          1 0.0000041 0.0000041   0.127    0.723    
W:L        1 0.0000107 0.0000107   0.336    0.565    
Residuals 45 0.0014317 0.0000318                     
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

> 
